VERSITA 


Chemical Papers 68 (12) 1725-1738 (2014) 
DOI: 10.2478/sll696-014-0622-y 


ORIGINAL PAPER 


Kinetics of thermal degradation of wood biomass 


Ivan Hrablay, Ludovit Jelemensky* 


Department of Chemical and Biochemical Engineering, Institute of Chemical and Environmental Engineering, 
Faculty of Chemical and Food Technology, Slovak University of Technology, Radlinskeho 9, 812 37, Bratislava, Slovakia 


Received 30 January 2014; Revised 9 June 2014; Accepted 10 July 2014 

Dedicated to the memory of professor Elemir Kossaczky 

Pyrolysis kinetics of a hardwood representative, beech (Fagus sylvatica ), was investigated by 
two different kinetic approaches: model-free isoconversional method and model-fitting method. The 
model-free isoconversional method was used for the determination of apparent kinetic parameters, 
i.e. the activation energy and pre-exponential factor. The model fitting method was used for the 
optimization of kinetic parameters of the reaction pathways of three selected reaction mechanisms: 
one-step, two-step, and three-step one. In both approaches, thermo-gravimetric data were used 
at five heating rates: 2°C min -1 , 5°C min -1 , 10 °C min -1 , 15°C min -1 and 20 °C min -1 . As the 
most suitable mechanism, the three-step mechanism containing the intermediate degradation step 
was chosen. This selection was supported by experimental results from the 13 C NMR analysis of 
solid residues prepared at the key temperatures within the range of 230-500 °C. The progress of 
mass fraction values of each component in this mechanism was simulated. Conclusions from the 
simulation were confronted with experimental results from the 13 C NMR. 
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Introduction 

In recent years, the need for renewable energy 
sources, especially biomass energy, has significantly 
increased due to the continual depletion of conven¬ 
tional fossil fuels. Another aspect is the air pollu¬ 
tion which is much lower considering the utilization of 
biomass energy in contrast to fossil fuels. There is also 
a growing pressure of government policies to achieve 
better environmental aspects of power generation pro¬ 
cesses. 

Wood and agricultural residues have a specific im¬ 
portance because of their wide distribution and easy 
access at relatively low costs. Wood is even more im¬ 
portant because of its higher density (higher energy 
content per volume), lower ash residue and very low 
nitrogen content. In principle, the transformation of 
energy from wood can be achieved by thermo-chemical 
processes: direct combustion, pyrolysis and gasifica¬ 


tion. However, the efficiency of energy transformation 
in direct combustion is low compared to that achieved 
by pyrolysis and gasification. Moreover, pyrolysis and 
gasification provide all three phases of products while 
each one of them can be maximized or minimized by 
process control. After applying certain cleaning pro¬ 
cedures, gaseous and liquid products from pyrolysis 
and gasification can be used in gas turbines or diesel 
engines. In addition, quality fuels such as methanol, 
hydrogen or fine chemicals can also be produced. The 
solid phase product represented by charcoal can also 
be used as quality fuel. 

Kinetic analysis 

Pyrolysis is the first step in any gasification or com¬ 
bustion process and thus understanding of pyrolysis 
kinetics is essential. Pyrolysis of solid materials, such 
as biomass, is considered as a heterogeneous chemical 
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reaction. More than thirty years ago, when the solid- 
state kinetic theory was still inconsistent, the study 
of cellulose pyrolysis mechanisms, the most abundant 
organic compound in nature and a representative of 
ligno-cellulosic materials, was only in its initial stage. 
Since then, the pyrolysis reaction scheme has evolved 
substantially and the reaction schemes for the whole 
ligno-cellulosic substrate have been improved (White 
et ah, 2011). However, it is well known that there are 
still many questions concerning the mechanism of py¬ 
rolysis because of the complexity of reactions and the 
variations in the physical and chemical properties of 
biomass. There is a great amount of factors affecting 
the pyrolysis rate, yields, composition and properties 
of the products. Temperature, pressure and heating 
rate are the main operating parameters. In addition, 
biomass properties such as chemical composition, in¬ 
organic salts content, particle size and shape, density, 
moisture content, etc. also play an important role in 
the pyrolysis process (Di Blasi, 2008). 

In general, kinetic studies involve process rates 
measurement and parameterization. The rate of the 
pyrolysis process can be parameterized to the form 
traditionally used for the biomass decomposition ki¬ 
netics (Eq. 1): 

^ = k(T)f(a) (1) 

where a is the extent of conversion, t is the time, k(T ) 
is the rate constant dependence on the temperature 
(T), and f(a) is the conversion function representing 
a reaction model and depending on the reaction con¬ 
trolling mechanism. 

The rate constant is typically expressed using the 
fundamental Arrhenius rate law (Eq. 2): 

k(T) = A exp ( 2 ) 

where A is the pre-exponential (frequency) factor, E 
is the activation energy of the reaction, and R is the 
universal gas constant. 

The extent of conversion, a , can be defined as the 
mass fraction of decomposed material (Eq. 3): 


where w is the actual mass of decomposing material, 
wo is the initial mass, and Wf is the final mass. 

The rate of the process in form of Eq. (1) describes 
the single-step process. If the process consists of mul¬ 
tiple steps, the process rate equation is more complex. 
For example, the process rate equation involving two 
parallel steps is expressed in form of Eq. (4): 

^ + k 2 {T)f 2 (a 2 ) (4) 

where oq and a 2 are the individual extents of conver¬ 
sion for a given reaction step, whereby their summa¬ 
tion gives an overall value, a. 


Kinetic analysis methods 

Principally, there are two types of methods used 
in kinetic analysis: the so called model-free and the 
model-fitting method. 

The model-free approach implies that the determi¬ 
nation of kinetic parameters can be done without as¬ 
suming a reaction model. Isoconversional methods are 
frequently used model-free methods (e.g. Vyazovkin & 
Dollimore, 1996; Flynn, 1997; Budrugeac, 2002; Burn¬ 
ham & Dinh, 2007; Sanchez-Jimenez et ah, 2010). Ac¬ 
cording to the isoconversional principle, the reaction 
rate at a constant extent of conversion is a function 
of temperature only and the resulting kinetic parame¬ 
ters are dependent on the conversion. The most com¬ 
monly used differential isoconversional method is the 
Friedman method (Friedman, 1964). For linear non- 
isothermal heating program, the equation derived by 
combining Eqs. (1) and (2) and the following lineariza¬ 
tion are commonly used (Eq. 5): 

= <5) 

where index j represents the given heating program, 
and index a represents a variable at the constant ex¬ 
tent of conversion. Input data consist of the extent 
of conversion at the given temperature while output 
data are represented by the pre-exponential factor and 
the activation energy at the given extent of conversion 
computed from the intercept and the slope of the so 
called isoconversional lines, respectively. It is neces¬ 
sary to choose the reaction model for the computation 
of the pre-exponential factor only. The most common 
form of the reaction model is: (1 - <a) n , where n is the 
reaction order. 

Model fitting methods use a particular reaction 
model to reach the best fit of experimental data re¬ 
sulting in the determination of kinetic parameters. 
Each model-fitting procedure implies minimization of 
the difference between the experimental and the calcu¬ 
lated reaction rate values. A wide variety of reaction 
models is listed in the project of the ICTAC kinetics 
committee (Vyazovkin et ah, 2011) together with ap¬ 
propriate model selection recommendations. Accord¬ 
ing to an ICTAC kinetic study (Brown et ah, 2000), 
reliability of the model-fitting methods is comparable 
with that of the model-free isoconversional methods if 
the fitting procedure is used simultaneously for multi¬ 
heating rates. A great advantage of the model-fitting 
methods, in contrast to the isoconversional ones, is 
the ability to treat multi-step degradation processes 
evaluating the respective kinetic constants of complex 
kinetics. 

The model fitting methods can be further divided 
into linear and non-linear. Non-linear methods are 
considered as more suitable than the linear ones be¬ 
cause the linearization causes distortion of kinetic pa- 
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rameters for the most important sequence of the pro¬ 
cess. Furthermore, the non-linear methods enable sim¬ 
ple optimization of process rates, conversions, or both 
simultaneously, and they are able to treat any set of 
differential rate equations using numerical integration 
(Vyazovkin et ah, 2011). 

Both single-step and multi-step processes can be 
fitted using the non-linear regression methods. An im¬ 
portant step of the model-fitting methods is the se¬ 
lection of an appropriate multi-step mechanism. The 
selection of a formal mechanism, such as parallel, con¬ 
secutive, reversible, or a combination of them, can 
be done based on the knowledge of possible mech¬ 
anisms of the investigated process. Another impor¬ 
tant step of model-fitting methods is the determina¬ 
tion of the number of individual steps present in a 
multi-step mechanism. Some information can be ob¬ 
tained from the experimental data; the reaction pro¬ 
file should be checked for the presence of shoulders 
and one or more peaks. The presence of a shoulder 
or a peak means that there is at least one reaction 
step in the process. However, the knowledge of poten¬ 
tial mechanisms helps significantly when determining 
the most adequate one. It has to be noted that there 
is a limitation of the number of reaction steps used 
due to computational complications and experimental 
data precision requirements. Computational compli¬ 
cations result from the interaction of too many E and 
A parameters so that a global minimum is difficult to 
find. It is recommended to add another reaction step 
only if this results in a significant decrease of the dif¬ 
ference between the experimental and the calculated 
data (Opfermann, 2000). In conclusion, for reliable 
kinetic parameters determination, non-linear model¬ 
fitting methods for multi-heating rates should be used. 
Treating the kinetics of individual components provid¬ 
ing, in summation, the kinetics of the whole reactant 
matter can be considered as a separate approach in 
the non-linear model-fitting methods (e.g. Gasparovic 
et ah, 2010; Yoon et ah, 2012; Haydary & Susa, 2013). 

Recently, Wu et al. (2013) published a simplified 
form of a non-linear model-fitting method to simulate 
the pyrolysis process of the cellulose model compound, 
glyceraldehyde. The authors proposed a two-step con¬ 
secutive kinetic model including an intermediate. The 
reaction path proposed includes direct fragmentation 
of glyceraldehyde to light volatiles and polymerization 
of glyceraldehyde to the intermediate (heavy volatiles) 
in the solid/liquid phase. In the second step, the in¬ 
termediate decomposes to char and volatiles. Each re¬ 
action step is described by its own rate equation. Fi¬ 
nally, the mass balance was formulated, balancing the 
actual mass of glyceraldehyde with that of volatile 
fractions released. The objective function (OF) was 
defined as the sum of square differences between the 
experimental and the calculated values of mass loss. 
Then, kinetic parameters for each pathway were de¬ 
termined and the mass fractions of each component in 


the mechanism were computed in dependence on the 
temperature. 

Pyrolysis mechanisms 

Biomass pyrolysis involves numerous complex re¬ 
actions resulting in a large number of intermediates 
and products; therefore, proposing an exact reaction 
mechanism and modeling the kinetic behavior is ex¬ 
tremely difficult. However, it can be done based on 
visible kinetics which motivated researchers to pro¬ 
pose a more precise mechanism employing available 
experimental techniques. However, even nowadays it 
is difficult to formulate a precise mechanism taking 
into account various feedstock and experimental con¬ 
ditions. All mechanisms can be described by a sim¬ 
plified general approach when biomass feedstock is 
transformed into gas/volatiles, tar, and char as the 
end products (Prakash & Karunanithi, 2008). Litera¬ 
ture presents various decomposition mechanisms ap¬ 
plicable to biomass pyrolysis. The general depiction 
of a biomass pyrolysis scheme includes moisture evap¬ 
oration followed by degradation of less stable poly¬ 
mers and, finally, decomposition of the most stable 
components which takes place while volatiles are re¬ 
leased (Arseneau, 1971; Fisher et al., 2002). Solid char 
is formed in the temperature range of 200-400 °C in 
the primary decomposition phase and it consequently 
reacts at temperatures over 400 °C, which is the sec¬ 
ondary pyrolysis stage (Fisher et al., 2002). 

Cellulose as the main component of biomass is of¬ 
ten used as its model compound because of its simpler 
structure. An incipient mechanism of cellulose pyrol¬ 
ysis, the two-step competitive mechanism describing 
cellulose transformation into volatiles in one step and 
char with gases in the second one was developed by 
Broido (Kilzer & Broido, 1965). According to the ex¬ 
perimental observations, Antal and co-workers (Antal 
& Varhegyi, 1995; Varhegyi et al., 1997; Antal et al., 
1998) proposed to use the one-step reaction mecha¬ 
nism of the first-order reaction equation to describe 
cellulose pyrolysis. Ultimately it was concluded (Varh¬ 
egyi et al., 1994) that complicated models including 
more than one reaction step are not necessary to simu¬ 
late the mass loss behavior of cellulose, justifying this 
with the rate limiting step, the depolymerization of 
cellulose. Depolymerization can be satisfactorily de¬ 
scribed by the one-step mechanism of the first-order 
reaction equation with high activation energy (Varh¬ 
egyi et al., 1994). However, this simplification seems 
too rough as the composition of light gases in pyrolytic 
mixtures (Banyasz et al., 2001a, 2001b; Li et al., 2001) 
proves that cellulose pyrolysis includes at least two 
reaction steps. These steps have to be competitive 
as their yields distinctly depend on the heating rate. 
Then the two-step competitive mechanism seems to be 
a rational and sufficient approximation of the pyroly¬ 
sis process. However, the hypothesis of two competing 
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reaction steps is not novel and it has been continu¬ 
ously modified and clarified (Mamleev et ah, 2007). 
Maybe the most important contribution has been done 
by inserting the so called “active cellulose” or “in¬ 
termediate” into the reaction scheme (e.g. Broido & 
Weinstein, 1972; Bradburry et ah, 1979; Shafizadeh, 
1982; Koufopanos et ah, 1989; Lede, 2012). This in¬ 
termediate is thought to be produced by an initial 
reaction step not associated with any mass loss while 
the degree of polymerization decreases. There is still 
no consensus on the existence of this intermediate in 
literature. Transformation of cellulose to active cellu¬ 
lose is reported as either very fast or not present at 
all (Varhegyi et ah, 1994; Antal & Varhegyi, 1995). In 
the three-step mechanism (Shafizadeh & Chin, 1977; 
Thurner & Mann, 1981) that assumes biomass decom¬ 
position via three parallel reactions, the intermediate 
is simply omitted. This could be an adequate overall 
representation, which is however not correct on the 
molecular level (Hoekstra et ah, 2012). 

Using modern analytical devices to analyze cel¬ 
lulose transformation reactions and their pyrolysis 
products seems to be very helpful in precise mecha¬ 
nism proposition. The first and most common prac¬ 
tice is to investigate pyrolysis decomposition kinet¬ 
ics by the thermo-gravimetric (TG) analysis, where 
a sample can be subjected to a wide range of con¬ 
trolled temperature programs. TG analysis in com¬ 
bination with other techniques of thermal analysis is 
an appropriate method to reach a mechanistic conclu¬ 
sion of the solid-state decomposition process (How¬ 
ell, 2006). Due to the high complexity of this process, 
also additional techniques have to be used to investi¬ 
gate structural and compositional changes of a given 
material. These techniques include the Fourier trans¬ 
form infrared spectroscopy (FTIR), gas chromatogra¬ 
phy/mass spectrometry (GC/MS), nuclear magnetic 
resonance (NMR), and others. The 13 C NMR tech¬ 
nique is more often used and cited in literature as a 
promising method providing important results regard¬ 
ing biomass transformation into char (e.g. Brewer et 
ah, 2009; Melkior et ah, 2012; Zheng et ah, 2013). 
These results enable the most appropriate degrada¬ 
tion mechanism proposition, and together with the 
data obtained from the thermal analysis, they serve 
as an input for the kinetic analysis providing ki¬ 
netic parameters of the degradation process. Shen and 
Gu (2010) investigated cellulose pyrolysis in detail 
using the analytical techniques Py-GS-MS and TG- 
FTIR proposing speculative chemical pathways of the 
main products formation. However, conclusions on the 
chemical pathways are very speculative due to the 
complex chemical structure of the pyrolysis products. 
Some efforts to better understand the chemical path¬ 
ways have been done using much simpler model com¬ 
pounds such as glucose, levoglucosan, and glycerine 
(e.g. Paine et ah, 2007; Hosoya et ah, 2008; Geng et 
ah, 2011). These compounds are formed as interme¬ 


diates or products during cellulose pyrolysis. Most re¬ 
cently, Wu et al. (2013) used glyceraldehydes, one of 
the products of cellulose pyrolysis, as a model com¬ 
pound as its chemical structure is more similar to that 
of glucose/cellulose than that of glycerine. This ap¬ 
proach of simplified model compounds is widely used 
although it should be kept in mind that these model 
compounds form only small percentages of the whole 
palette of pyrolysis products. 

In this paper, the model-free isoconversional me¬ 
thod was used to acquire an initial insight into the 
process complexity of wood biomass pyrolysis and to 
determine apparent kinetic parameters. Consequently, 
the non-linear model-fitting method was used with the 
initial estimates of kinetic parameters obtained by the 
isoconversional method to fit the experimental data 
using various reaction mechanisms. Discrimination of 
these mechanisms was done based on the experimental 
results from the 13 C NMR analysis of solid residues 
prepared at key temperatures. 

Experimental 

Experiments were conducted using powder sam¬ 
ples of hardwood, beech (Fagus sylvatica ), prepared by 
milling in a ball mill. No thermal pretreatment was ap¬ 
plied. The mass of samples was approximately 20 mg 
each. Experiments were conducted in a simultaneous 
thermo-gravimetric analyzer STA 409 PC Luxx (Net- 
zsch). Reaction environment was purged with argon at 
60 mL min -1 . For kinetic parameters determination, 
the experiments were conducted starting from ambi¬ 
ent temperature up to 600 °C using five heating rates: 
2°C min -1 , 5°C min -1 , 10 °C min -1 , 15 °C min -1 and 
20 °C min -1 . For the solid residues preparation, exper¬ 
iments started from ambient temperature up to the 
temperatures listed in Table 1. The temperature pro¬ 
gram consisted of a dynamic sequence with the heating 
rate of 5°C min -1 and an isothermal sequence applied 
for 20 min. Temperatures of samples preparation were 
selected according to the TG curve characterizing its 
important points. In this way, samples prepared at the 
selected temperatures present the solid residue charac¬ 
terizing structural changes in wood matter during the 
pyrolysis up to a certain extent. Table 1 includes also 
the mass fractions of C and H obtained by elemental 
analysis of solid residues used for the NMR analysis. 
Values of the mass fractions of C, H, and that of the 
C/H ratio only are published. 

NMR analysis 

Solid residues prepared were analyzed using a Var- 
ian VNMRS 600 MHz (Oxford Instruments) equipped 
with a superconducting magnet of the induction of 
14.1 T, and a 3.2 mm probe (Narrow Bore Triple Res¬ 
onance HXY MAS). The samples were placed in thick- 
wall rotors with the rotation of 18 kHz. The correct 
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Table 1 . Temperatures of samples preparation with the corresponding mass loss and mass fractions of C and H obtained by 
elemental analysis 


Final temperature 

Mass loss 


Mass fractions and its ratio a 


C 

H 

C/H 

°C 

% 


mass % 

- 

230 

11.8 

47.9 (0.0) 

6.0 (0.7) 

7.9 

250 

18.2 

50.1 (0.1) 

7.0 (0.7) 

7.2 

270 

28.9 

50.7 (0.1) 

5.9 (0.0) 

8.6 

290 

40.1 

54.7 (0.2) 

5.7 (0.0) 

9.6 

300 

44.5 

- 

- 

- 

310 

58.7 

61.0 (0.1) 

5.2 (0.0) 

11.8 

320 

63.1 

- 

- 

- 

330 

64.8 

69.2 (0.1) 

4.5 (0.0) 

15.4 

340 

66.2 

- 

- 

- 

350 

72.0 

70.9 (0.3) 

4.1 (0.1) 

17.5 

410 

77.5 

72.8 (0.1) 

3.3 (0.1) 

21.9 

500 

77.5 

74.3 (0.1) 

3.1 (0.1) 

23.6 


a) Values for C and H are mean values, values in parentheses are the absolute deviations computed from repeated analyses. 


a b 




Fig. 1 . TG (a) and DTG (b) curves of beech pyrolysis at five heating rates and the argon flow rate of 60 mL min 1 . Legend: 
2°C min” 1 (□), 5°C min" 1 (O), 10°C min” 1 (A), 15°C min” 1 (V), 20°C min" 1 ( 0 ). 


spectral scale was determined using the CH 2 group of 
adamantine as the reference. Impulse sequence con¬ 
sisted of: relaxation -5 s, spin echo, and acquisition - 
19.3 ms. 

Results and discussion 

TG curves obtained at the heating rates of 
2°C min” 1 , 5 °C min -1 , 10 °C min” 1 , 15 °C min” 1 
and 20 °C min -1 , respectively, and the first derivative 
of the TG curves, the DTG curves, are depicted in 
Fig. 1. As it can be seen from the DTG curves, ther¬ 
mal degradation shows three peaks. The first one, up 
to about 160 °C, corresponds to sample drying with 
the mass loss of about 5 mass %. The second one, in 
the temperature interval of 250-320^0, represents the 
degradation mainly of hemicellulose. The third one, 
in the temperature interval of 320-400 °C, represents 
the degradation mainly of cellulose. The temperature 


interval of 250-400 °C represents active pyrolysis with 
the mass loss of about 70 mass %. Over 400 °C, passive 
pyrolysis takes place. However, lignin degradation oc¬ 
curs in the broad temperature interval of 200-500 T! 
and cannot be distinguished from hemicellulose and 
cellulose degradation. The peak complexity, more ob¬ 
vious at higher heating rates, is due to the relatively 
high content of hemicellulose degrading at lower tem¬ 
peratures as cellulose. Considering the elemental com¬ 
position, it is obvious that the fraction of C increases 
with the increasing temperature at the expense of the 
decreasing fraction of H. This trend confirms that the 
solid residue is more charred due to the increasing 
temperature. 

Results from the 13 C NMR analysis are shown 
in Fig. 2. Chemical shifts of cellulose, hemicellulose, 
and lignin groups were assigned according to Table 2 
based on many studies published over the last decades 
(Bardet et ah, 1997; Gil & Neto, 1999; Maunu, 2002). 
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Table 2. Chemical shifts of carbon identified in the spectra of solid residues prepared at the temperatures in Table 1 


Chemical shift, (5 

Carbon specification 

Labeling 

21 

Hemicellulose CH 3 COO— 

HI 

56 

Lignin —OCH 3 

LI 

62 

Cellulose C -6 (amorphous) 

Cl 

65 

Cellulose C -6 (crystalline) 

C 2 

72 

Cellulose C-2/3/5 

C3 

75 

Cellulose C-2/3/5 

C4 

84 

Cellulose C-4 (amorphous) 

C5 

89 

Cellulose C-4 (crystalline) 

C 6 

105 

Cellulose C-l 

C7 

121 

Lignin G -6 

L2 

133-134 

Lignin S-l (ne), S-4 (ne) Lignin G-l (ne) 

L3 

138-139 

Lignin S-l (e), S-4 (e) Lignin G-l (e) 

L4 

147-148 

Lignin S-3 (ne), S-5 (ne) Lignin G-3 (ne, e), G-4 (ne, e) 

L5 

153 

Lignin S-3 (e), S-5 (e) 

L 6 

172 

Hemicellulose -COO-R Hemicellulose CH 3 -COO- 

H2 


S: carbon in syringyls, G: carbon in guaiacyls, ne: in non-etherified arylglycerol /3-aryl ethers, e: in etherified arylglycerol /3-aryl 
ethers. 



Fig. 3. Carbon position labeling in the main wood compo¬ 
nents: cellulose, hemicellulose (acetyl groups are not 
shown due to the different hydroxyl group substitution) 
(a), lignin monomer (guaiacyl - methoxyl group at po¬ 
sition 3, syringyl - methoxyl groups at positions 3 and 

5) (b). 


The position of carbon was labeled according to Fig. 3. 

Comparing the NMR spectra obtained at differ¬ 
ent temperatures, two distinct temperature sections 
can be identified: under and over 300 °C. Spectra in 
the first section seem to be identical; however, de¬ 
tailed inspection revealed some changes. The hemi¬ 
cellulose acetyl methyl (HI) signal continuously de¬ 
creases with the increasing temperature, which means 
that hemicellulose is being degraded. The hemicel¬ 
lulose acetyl carboxyl (H2) signal decrease also cor¬ 
responds with this statement. The highest tempera¬ 
ture with detectable hemicellulose signal is 320 °C. The 
cellulose signals in the spectra are the most obvious 
ones. Cellulose C-6 (amorphous) (Cl) and cellulose C- 


4 (amorphous) (C5) signals representing amorphous 
cellulose are both more evident at 300 °C while cel¬ 
lulose C-6 (crystalline) (C2) and cellulose C-4 (crys¬ 
talline) (C6) signals representing crystalline cellulose 
start to decrease between 270-290 °C. This evidently 
shows the transformation of the crystalline form into 
an amorphous one. The cellulose C-l (C7) signal rep¬ 
resenting cellulose as a whole remains unchanged up to 
300 °C. Cellulose is the most degraded in the temper¬ 
ature range of 300-310 °C. Only small residues remain 
at 320 °C. The second most obvious signals in the spec¬ 
tra are those of lignins. The lignin methoxyl (LI) sig¬ 
nal is quite persistent and becomes more distinguished 
as the cellulose degrades at 320 °C. With further tem¬ 
perature increase, this signal continuously decreases, 
which represents the process of lignin demethoxyla- 
tion. The lignin G-6 (L2) signal increases very slightly 
together with the aromatic lignin signals in the chem¬ 
ical shift range of 115-160 ppm, and individually from 
310 °C. As it represents the guaiacyl structure, it can 
be concluded that its increasing signal is caused by 
lignin demethoxylation due to its transformation from 
the syringyl structure. This fact corresponds with the 
previous statement. The lignin S-l (ne), S-4 (ne) and 
lignin G-l (ne) (L3) signal is also of importance, not 
because it represents both the syringyl and the gua¬ 
iacyl structure but because both structures are in 
their non-etherified form. An increase of this signal 
then necessarily means cleavage of the (3-0 -4 bonds 
between lignin dimers. This behavior can be clearly 
seen at above 290°C. The lignin S-l (e), S-4 (e) and 
lignin G-l (e) (L4) signal representing the etherified 
forms is then expected to decrease. Although this sig¬ 
nal is quite low, its decrease can be observed between 
270-290 °C. Behavior of the last two lignin signals fits 
perfectly into the already drawn mosaic. The lignin 
S-3 (ne), S-5 (ne) and lignin G-3 (ne, e), G-4 (ne, 
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e) signal (L5) representing both the etherified and 
the non-etherified forms increases very distinctly from 
the lowest temperature up to 300 °C, then it remains 
constant up to 320 °C when it rapidly decreases. At 
the same time, the lignin S-3 (e), S-5 (e) signal (L6) 
representing the etherified form, seems to be merged 
with the increasing L5 signal. However, this behav¬ 
ior can represent either the (5-0 -4 bonds cleavage or 
the lignin demethoxylation into guaiacyl. Two follow¬ 
ing facts support the lignin demethoxylation pathway: 
1) (5-0 -4 bonds cleavage is observed at temperatures 
not lower than 270-290 °C, and 2) lignin G-6 (L2) sig¬ 
nal slightly increases already from the lowest tempera¬ 
ture. However, starting from the temperature range of 
270-290 °C, also the (5-0 -4 bonds cleavage takes place, 
which results in a much intensive L5 increase. The 
best readable lignin structure can be seen at about 
320 °C. Starting from 320°C, not only the LI signal 
decrease due to lignin demethoxylation contributes to 
the course of lignin degradation but also the L5 signal 
rapidly decreases. As at this temperature, the L5 sig¬ 
nal is supposed to consist of pure guaiacyl structure 
in etherified or non-etherified form, its decrease means 
either its depolymerization into monomer units or the 
degradation of the guaiacyl monomer unit itself. At 
the last two temperatures, 410°C and 500 °C, signals 
in the chemical shift range of 110-140 ppm are abso¬ 
lutely dominant making the original lignin structure 
to disappear. At this point, all that remains is the 
charred residue containing condensed aromatic struc¬ 
tures originating from all wood components. 

In conclusion, hemicellulose acetyl groups release 
is observed at above 230 °C, the last hemicellulose 
residues can be found at 320 °C. Regarding cellulose, 
its transformation from the crystalline form into an 
amorphous one starts in the temperature range of 
270-290 °C and continues up to 310 °C. The highest 
cellulose degradation occurs between 300-310 °C, only 
small residues remain at 320 °C. Regarding lignin, at 
lower temperatures, of up to 270°C, the process of 
lignin demethoxylation starts preferentially. Cleavage 
of the (5- O -4 bonds between lignin dimers starts to oc¬ 
cur between 270-290 °C. With the increasing temper¬ 
ature, the combination of both processes takes place. 
The final process includes the formation of condensed 
aromatic structures and the solid residue charring. 

Model-free isoconversional method 

The model-free differential isoconversional method 
was used to provide the apparent kinetic parameters. 
Input data are in form of TG curves at five differ¬ 
ent heating rates. The temperature range was cut to 
160-500 °C excluding sample drying and the final tem¬ 
peratures where negligible mass loss was observed. In 
addition, the conversion range was cut to 10-75 % 
(8-57 % of mass loss) due to the kinetic parameters 
uncertainty intrinsic to this method. The conversion 



Fig. 4. Isoconversional lines for predefined conversions in the 
range of 10-75 % for five heating rates. 



Fig. 5. Dependence of the activation energy on the conversion 
computed from the isoconversional lines. 



Fig. 6. Dependence of the pre-exponential factor on the con¬ 
version computed from the isoconversional lines. 


range was divided using a step of 0.005. Eq. (5) is visu¬ 
ally interpreted by the isoconversional lines presented 
in Fig. 4. The slope and intercept of these lines en¬ 
able the computation of the activation energy (Fig. 5) 
and the pre-exponential factor (Fig. 6), respectively. 
The activation energy and pre-exponential factor at 
selected conversions with their mean values are listed 
in Table 3. 

As it can be seen from Figs. 5 and 6, and from 
Table 3, activation energy changes in the range of 
190-240 kJ mol -1 (mean value of 226.2 kJ mol -1 ) 
and the pre-exponential factor in the range of 10 15 - 
10 19 s -1 (mean value of 1.90 x 10 18 s -1 ). Such param- 









I. Hrablay, L. Jelemensky/Chemical Papers 68 (12) 1725-1738 (2014) 


1733 


Table 3. 

Kinetic parameters 
mean values 

at selected conversions and their 

a 

E 


A 


- 

kJ mol -1 



L 

0.1 

190.4 

1.09 

X 

10 15 

0.2 

211.2 

4.98 

X 

10 16 

0.3 

237.5 

4.95 

X 

10 18 

0.4 

242.1 

4.87 

X 

10 18 

0.5 

232.1 

3.34 

X 

10 17 

0.6 

226.7 

7.94 

X 

10 16 

0.7 

224.2 

3.93 

X 

10 16 

0.75 

225.2 

4.28 

X 

10 16 

Mean® 

226.2 

1.90 

X 

10 18 


a) Mean values were computed using each point from the con¬ 
version range. 


A 

B 

C 


wood —~—► vvolatiles + char 


woodl - 

, 

► v,vo(atilesl 


wood2 - 


► v 2 volatiles2 + char 


woodl - 


* Vj volatilesl 


wood2 - 

an 

* intermediate— 

-► v ? volatiles2 + char 


Fig. 7. Selected mechanisms of wood pyrolysis for the model¬ 
fitting procedure. 


V 2 portion of volatiles2. Mechanism C is the three- 
step mechanism describing wood transformation into 
an intermediate in the second parallel step. The third 
step describes the intermediate’s transformation into 
char. This mechanism is the extension of the previ¬ 
ous one with an additional component, the intermedi¬ 
ate. The original extension of the Kilzer-Broido mech¬ 
anism was presented by Shafizadeh and co-workers 
(Bradbury et ah, 1979); where the intermediate was 
produced directly from cellulose without any volatiles 
being released. Mechanism C in the form seen in Fig. 7 
has been presented only recently (Wu et ah, 2013). 

As the portions of released volatiles are unknown, 
the only experimental variable that can be compared 
against the computed one is the mass loss, i.e. the TG 
curve. The mass of the sample according to reaction 
mechanisms A, B, and C can be defined by Eqs. (6), 
(7) and (8), respectively. Generally, the sample mass 
equals to the original sample mass minus the volatiles 
with the corresponding portions, where aq, a \2 and 
c \2 are conversions of the given reaction step. 

^sample = ^original ^volatiles = 

= m-original UHUT-original (b) 

^sample — ^original ^volatilesl ^volatiles2 — 

= Unoriginal Ui O'i Unoriginal ^2 0^2 Unoriginal (7) 

^sample : ^original ^volatilesl ^volatiles2 = 

= Unoriginal ^lOqUn 0 riginal t>2 0^2 ^intermediate (8) 


eter changes indicate a complex process. Both, activa¬ 
tion energy and pre-exponential factor dependences 
exhibit a very similar trend. With the increasing con¬ 
version (temperature), they increase rapidly reaching 
the maximum in the conversion range of 30-40 % (23- 
31 % of mass loss), which corresponds to the temper¬ 
ature range of about 300-320 °C. With a further con¬ 
version increase, both parameters gradually decrease. 
Their mean values were used as the initial values for 
the model-fitting methods. 


The reaction model used was in the form of 
(1 - a) n and the kinetic constants in the Arrhenius 
form. The rate equation according to reaction mecha¬ 
nism A has the form of Eq. (9), while reaction mech¬ 
anisms B can be described by Eqs. (10) and (11), and 
reaction mechanism C by Eqs. (10), (12) and (13). 

^ =Aexp (^) <i ~ a>n <9> 

t L = Aiexp (i§-) <1 - a ‘ r (10) 


Model-fitting method 

For the model-fitting procedure, the three mecha¬ 
nisms presented in Fig. 7 were selected. 

Mechanism A is the one-step mechanism describ¬ 
ing wood transformation into char while volatiles are 
released. Symbol v represents the final mass fraction 
of volatiles. This mechanism was selected as the sim¬ 
plest one used in the initial stages of modeling (e.g. 
Miyanami et ah, 1977; Wichman & Atreya, 1987; 
Antal & Varhegyi, 1995). Mechanism B is the two- 
step mechanism describing wood transformation in 
two parallel steps. The first step produces v\ portion 
of volatilesl, the second one produces char releasing 


^r = - 4jexp (i!0 (1 -“ j) ” (n) 

^ 3 r & ' 4l!exp (^ i ) (1 ““ i2) “ (12) 

= A 2 exp (“12 - “2)" 2 (13) 

Parameters to be optimized are: A, E : n and v, in¬ 
dexed according to the mechanism used. Integration of 
the rate equations provides the respective conversions. 
Computed sample mass has to be recalculated to the 
mass loss derivative entering the objective function 
(OF). The OF to be minimized is defined as a sum of 
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Table 4. Optimized parameters for three selected mechanisms with the corresponding OF values 


Mechanism A 


A/s- 1 

6.00 x 10 7 

E /(kj mol- 1 ) 

111.1 

n 

1.85 

V 

0.76 

OF 

6.68 x 10 -5 


Mechanism B 

Ai/s- 1 

7.48 x 10 15 

A 2 /s~ 1 

2.72 x 10 19 

£i/(kJ mol- 1 ) 

190.7 

E 2 /{kJ mol- 1 ) 

250.0 

m 

3.00 

n 2 

1.80 

Vl 

0.40 

V2 

0.31 

OF 


2.80 x 10 -5 



Mechanism C 


An/s 1 

7.60 x 10 11 

Ai 2 /s 1 

8.15 x 10 19 

A 2 /s 1 

2.98 X 10 19 

Ei/(kJ mol -1 ) 

148.8 

Fi 2 /(kJ mol -1 ) 

210.9 

E 2 /(kJ mol -1 ) 

250.0 

nn 

1.45 


1.55 

n 2 

2.02 

Vl 

0.33 

ni2 

V2 

0.39 

OF 



2.83 x 10 -5 




a 



b 



c 



Fig. 8 . Experimental and calculated da/d t curves for mechanism A (a), B (b) and C (c) for heating rates of 2°C min -1 , 5°C min -1 , 
10°C min -1 , 15°C min -1 and 20°C min -1 (from above). Legend: selected mechanism (line), experimental data (dashed 
line). 
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square differences between the experimental and the 
calculated mass loss time derivatives considering all 
heating rates (Eq. 14): 



where p is the total number of points for the given 
heating rate /3. 

Optimization was performed in the Matlab compu¬ 
tational software. Integration of the rate equations was 
performed using the odel5s function. Minimization of 
the OF was performed using the fmincon function. To 
ensure that the global minimum of the OF is reached, 
the procedure of various initial estimates was used. 

The resulting optimized parameters are summa¬ 
rized in Table 4. As it can be seen, v in mechanism 
A, and the sum of v± and v 2 in mechanisms B and 
C approach the value of 0.76, which is the mean ex¬ 
perimental value of all five heating rates. Comparing 
the kinetic parameters obtained by the isoconversional 
method, mechanism A shows much lower values of 
these parameters. This means that mechanism A is 
inconvenient for the description of the process com¬ 
plexity. On the other hand, comparing these param¬ 
eters with those obtained by mechanisms B and C, 
similar values of both E and A can be seen at least for 
one reaction step; the values are very close to the sec¬ 
ond reaction step of mechanism B and C. The values 
of OF in mechanisms B and C are lower than one half 
of that in mechanism A, which confirms the relevance 
of at least the second reaction step. Introducing the 
intermediate step does not result in a lower value of 
OF in mechanism C; however, it is very similar to that 
in mechanism B. In this case, to justify the necessity 
of another reaction step, additional information from 
experiments is needed. 

Optimized parameters were used to solve Eq. (6), 
(7), and (8). Consequently, conversions and their 
derivatives were computed and compared with the ex¬ 
perimental data. Experimental and calculated da/dt 
curves for all heating rates are presented in Fig. 8. 

As it can be seen, mechanism A does not fit the 
experimental data at all. Mechanisms B and C fit the 
experimental data satisfactorily, while mechanism C 
fits higher heating rates better. 

Mass fractions of each component in mechanism C 
were determined in dependence on the temperature 
at the heating rate of 20°C min -1 (Fig. 9). For bet¬ 
ter comparison with the NMR results, gridlines corre¬ 
sponding to the temperatures from Table 1 are shown. 
This depiction has to be interpreted as mass fractions 
of each component at the given temperature instead 
of the time trend of these mass fractions. 

As it can be seen in Fig. 9, the mass fraction of the 
intermediate starts to increase followed by volatiles 1 
with the increasing temperature and the decreasing 



Fig. 9. Mass fractions of components in mechanism C (20 °C 
min -1 ) with gridlines corresponding to the tempera¬ 
tures from Table 1. Legend: experimental (•), sam¬ 
ple (□), volatiles 1 (O), volatiles2 (A), intermediate (V), 
char (0). 


sample mass. The intermediate is present in the tem¬ 
perature range of 200-400 °C; during its decomposi¬ 
tion, volatiles2 are released forming char. The final 
mass fraction of volatiles 1 is 0.33 and that of volatiles2 
is 0.39. 

Confrontation of simulation and NMR results 

As a result of the optimization procedure, mecha¬ 
nism B and C have shown approximately the same OF 
value. Mechanism C is slightly more convenient to fit 
the experimental data at higher heating rates. How¬ 
ever, discrimination between these two mechanisms 
has to be done using experimental results from the 
NMR analysis. 

According to the simulation results, the only 
change that can be seen below 230 °C, is the formation 
of the intermediate. In the NMR spectra, this temper¬ 
ature represents only the process of lignin demethoxy- 
lation. So it can be concluded that the presence of the 
intermediate is connected with the demethoxylation 
of the lignin structure. In the temperature interval of 
230-300 °C, changes mainly connected with the release 
of volatiles 1 should occur. Considering the release 
of volatiles only, this temperature interval is charac¬ 
terized by the release of hemicellulose acetyl groups 
and also by the process of lignin demethoxylation. At 
320 °C, the last residues of hemicellulose and cellulose 
can be found, while volatiles2 should be present. This 
could mean that the release of volatiles2 is not con¬ 
nected with the release of hemicellulose acetyl groups. 
On the other hand, at this temperature, the process 
of further lignin demethoxylation starts, which im¬ 
plies that the release of volatiles2 is connected only 
with the lignin demethoxylation. Along with the re¬ 
lease of volatiles2, char should be formed. In NMR, 
this can be accompanied with the cleavage of /3-0 -4 
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Fig. 10. Visual interpretation of simulation and NMR analysis 
results with additional information from the TG anal¬ 
ysis in dependence on temperature (c - crystalline, 
a - amorphous form of cellulose). 


bonds or the lignin monomer unit degradation result¬ 
ing in lignin-derived aromatic structures. The forma¬ 
tion of the main portion of char formed is accompanied 
with the formation of condensed aromatic structures. 
The intermediate degradation is finished before reach¬ 
ing the temperature of 400 °C and the presence of the 
intermediate should not be accompanied with the for¬ 
mation of condensed aromatic structures. 

According to the aforementioned facts, the inter¬ 
mediate represents a solid structure being neither a 
raw sample nor charred residue. This structure con¬ 
tinuously undergoes the transformations connected 
with the lignin structure degradation. The whole pro¬ 
cess described including the degradation intervals of 
biomass components and the position of the inter¬ 
mediate is visualized in dependence on temperature 
(Fig. 10). 


Conclusions 

In this paper, two different kinetic approaches 
to the determination of kinetic parameters of wood 
biomass pyrolysis were used. Results from the isocon- 
versional method confirm the complexity of the pro¬ 
cess. The so obtained mean value of the activation en¬ 
ergy is 226.2 kJ mol -1 and that of the pre-exponential 
factor is 1.90 x 10 18 s -1 . Optimization of kinetic pa¬ 
rameters according to three selected mechanisms was 
performed. Mechanism C was found to provide slightly 
better fit of the experimental data than mechanism B 
with comparable OF values. Although mechanism C 
provides a better fit of the peak complexity, its se¬ 
lection has to be confronted with experimental NMR 
results of solid residues prepared at the key tempera¬ 
tures of the TG curve. NMR results confirm the for¬ 
mation of an intermediate; its decomposition into char 


and volatiles2 does not affect the original raw sam¬ 
ple, nevertheless it is a continuously changing solid 
structure. Mass fractions of each component in mech¬ 
anism C were computed as the function of tempera¬ 
ture and confronted with the NMR results. It can be 
concluded that the presence of the intermediate is con¬ 
nected with the lignin structure degradation starting 
with the process of lignin demethoxylation. The for¬ 
mation of volatiles 1 is preferably connected with the 
release of hemicellulose acetyl groups and also with the 
process of lignin demethoxylation. On the other hand, 
the formation of volatiles2 is not connected with the 
release of hemicellulose acetyl groups but with the pro¬ 
cess of further lignin demethoxylation. The formation 
of char can be accompanied with the cleavage of /3- 0-4 
bonds or the lignin monomer unit degradation result¬ 
ing in lignin-derived aromatic structures. The forma¬ 
tion of the main portion of char is accompanied with 
the formation of condensed aromatic structures. The 
presence of an intermediate should not be accompa¬ 
nied with the formation of condensed aromatic struc¬ 
tures. The intermediate is a solid structure, neither a 
raw sample nor charred residue, which continuously 
undergoes transformations due to the lignin structure 
changes. 
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